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ABSTRACT 

The X-ray properties of a relaxed cluster of galaxies are determined primarily 
by its gravitational potential well and the entropy distribution of its intracluster 
gas. That entropy distribution reflects both the accretion history of the clus- 
ter and the feedback processes which limit the condensation of intracluster gas. 
Here we present Chandra observations of the core entropy profiles of nine classic 
"cooling-flow" clusters that appear relaxed and contain intracluster gas with a 
cooling time less than a Hubble time. We show that those entropy profiles are 
remarkably similar, despite the fact that the clusters range over a factor of three 
in temperature. They typically have an entropy level of ~ 130 keV cm^ at 100 kpc 
that declines to a plateau ~ lOkeVcm^ at < 10 kpc. Between these radii, the 
entropy profiles are oc r" with a ^ 1.0 — 1.3. The non-zero central entropy levels 
in these clusters correspond to a cooling time ~ 10^ yr, suggesting that episodic 
heating on this timescale maintains the central entropy profile in a quasi-steady 
state. 

Subject headings: catalogs - galaxies: clusters: general - X-rays: galaxies: clus- 
ters - cosmology: observations - methods: data analysis 

1. Introduction 

The global properties of a cluster of galaxies, such as its bolometric X-ray luminosity Lx 
and its mean temperature Tx , are determined primarily by the mass M^h within a suitably 
chosen virial radius. A cluster's temperature depends on mass because mass determines the 
depth of the cluster's potential well. Its X-ray luminosity depends on mass because mass 
determines both the total number of baryons in the cluster and the potential well confining 
those baryons. However, several secondary factors combine to produce a dispersion in both 
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Lx and Tx at a fixed Myir, and understanding the nature of that dispersion is crucial to 
doing precision cosmology with clusters. One of those factors is merger shocks, which can 
temporarily raise both the luminosity and best-fitting temperature of a cluster (e.g., Randall 
et al. 2002). A second is the shape of the potential well, because clusters whose potentials 
are more centrally concentrated tend to have higher central temperatures (e.g., Voit et al. 
2002). A third factor is the amount of intracluster gas with a cooling time less than the 
age of the universe. The presence of such gas leads to both a large peak in the central 
surface brightness of a cluster and a central temperature gradient that rises with radius. 
Consequently, clusters having larger amounts of gas with a short cooling time tend to have 
higher Lx and lower Tx at a given value of Myir (Allen & Fabian 1998; Fabian et al. 1994; 
Markevitch 1998). 

Such clusters have often been called "cooling-flow clusters," because the central gas was 
thought to condense and flow toward the center of the cluster as it radiated away its thermal 
energy (see Donahue & Voit 2004, for a recent review). Observations from Chandra and 
XMM-Newton now show that the central gas is not simply cooling to low temperatures and 
condensing in the manner originally envisioned (e.g., Peterson et al. 2001, 2003). Some form 
of feedback apparently prevents the central gas from condensing and forming stars, thereby 
truncating the high end of the galaxy luminosity function. The nature of that feedback is 
currently an active topic of both observational and theoretical research, focusing largely on 
the role of outflows from active galactic nuclei in cluster cores. 

This paper analyzes archival Chandra data on nine cooling-fiow clusters seeking clues 
to what keeps that gas from condensing and why clusters of a given mass have different 
amounts of gas with a short central cooling time. The tactic we take in our analysis is 
to focus on the entropy profiles of these clusters. We concentrate on entropy because it 
is a more fundamental property of the intracluster medium itself than either temperature 
or density alone. For example, the temperature of a cluster's gas primarily reflects the 
cluster's potential well depth; heating or cooling of the gas merely causes it to expand or 
contract in the potential well with only a modest change in temperature. The density of 
that gas depends on how much gravity can compress it in the cluster's potential well, and it 
is the specific entropy of the gas that determines its density at a given pressure. Thus, the 
observable X-ray properties of a relaxed cluster of galaxies depend almost entirely on two 
physical attributes: (1) the shape and depth of the cluster's dark-matter halo, and (2) the 
entropy distribution of the intracluster gas. (e.g., Voit et al. 2002). 

Intracluster entropy is also intimately related to the cooling and feedback processes that 
govern galaxy evolution and that may also play a role in limiting condensation in cluster 
cores. Theories and simulations of cluster formation which ignore these processes fail to 
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reproduce the observable properties of present-day clusters. If gravity alone were responsible 
for shaping the appearances of clusters and groups, then we would expect their properties 
to be nearly self-similar, with a luminosity-temperature relation like L oc T^. Furthermore, 
we would expect groups and clusters to have similar surface-brightness profiles, when scaled 
to the virial radius of the system. However, observations indicate that L oc (Edge 
& Stewart 1991; Markevitch 1998) and that the surface-brightness profiles of groups are 
shallower than those of clusters (Horner et al. 1999; Ponman et al. 1999). 

Many papers in the literature have attributed these deviations from self-similarity to 
an early episode of preheating by supernovac or active galaxies. Heat input that establishes 
a uniform minimum entropy level in the intergalactic medium breaks self-similarity because 
the extra entropy makes that gas harder to compress as it accretes into dark-matter haloes 
(e.g., Kaiser 1991; Evrard & Henry 1991). Thus, gas in the smaller potential wells of groups 
would be less centrally concentrated than gas in the larger potential wells of clusters, thereby 
reducing the himinosities of groups relative to clusters and flattening their surface-brightness 
profiles. Such a model predicts that the central entropy profile would fiatten into a core 
refiecting the minimum entropy level of the IGM. However, several analyses have suggested 
that the energy input required to explain the observed relations through global preheating 
is implausibly extreme. 

Voit & Bryan (2001) have argued that the entropy scale responsible for similarity break- 
ing is not a global property of the intergalactic medium but rather a requirement set by 
radiative cooling — the observed entropy at the core radii of groups and clusters turns out 
to be similar to the entropy level at which intracluster gas would cool within a Hubble 
time. Gas below this cooling threshold is therefore subject to cooling, condensation, and 
whatever feedback follows from that condensation, which may include both supernova and 
AGN activity. In this scenario, cooling and feedback conspire to deplete the amount of gas 
below the cooling threshold: some of the low-entropy gas condenses and feedback subse- 
quently raises the entropy of the remaining gas until both cooling and feedback shut down. 
Voit et al. (2002) have built upon this principle to create cluster models that account for 
many observable X-ray properties, including the Lx - Tx relation, the M^ir - Tx relation, 
the surface-brightness profiles of clusters, and their temperature gradients. However, those 
models do not explain why clusters differ in the amount of gas that still remains below the 
coohng threshold. 

Cluster entropy profiles have been presented by previous workers, with some of the 
earliest work coming from ROSAT and ASCA (Ponman et al. 1999; Lloyd-Davies et al. 
2000). They pointed out that if gravitational collapse was the only physical process affecting 
the intracluster gas, all entropy profiles should be self-similar. The temperature and density 
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profiles available from ROSAT and ASCA data were very limited, however. Results from 
XMM (Piffaretti et al. 2005) show that outside their core radii, clusters are similar, but the 
XMM data for inside the core radius of a cluster are limited by the spatial resolution of the 
telescope. 

In order to better understand the sub-threshold gas, the processes that keep some of it 
from condensing, and the effects of those processes on the global X-ray properties of clusters, 
we are developing a Chandra library of intracluster core entropy distributions. We selected 
the first batch of clusters for the library using criteria designed to ensure high-quality radial 
entropy profiles. The number of X-ray events had to be sufficient for extraction of high- 
quality spectra in multiple annular regions, and the clusters had to be relaxed enough to 
have an unambiguous central emission peak. We therefore compiled a master list of clusters 
available in the Chandra public archive as of December 1, 2003, when this project began, and 
cross-correlated it with the NASA/IPAC Extragalactic Database (NED) to obtain redshifts 
and positional data for the clusters. Wc then added any observations noted as cluster 
observations (those with sequence numbers beginning with '8') that did not contain NED 
clusters. After compiling this master list, we created Chandra images of all of the clusters 
and examined them to eliminate any clusters with insufficient counts for our study or with 
obvious substructure, such double peaks. We also eliminated observations badly affected by 
background flares. 

This selection process yielded a sample consisting entirely of nearby "cooling flow" clus- 
ters, most of which arc considered classic examples of that category. Table 1 lists that sample 
of nine clusters, including the cluster name, Chandra observation identification numbers, co- 
ordinates, maximum extraction radius, the redshift of the cluster used in spectral fits, its 
X-ray luminosity, and its mean X-ray temperature. These last two quantities are taken from 
ASCA observations analyzed by D. Horner in his PhD thesis.^ 

That is why this paper describing the first installment in our Chandra entropy library, 
which will ultimately contain clusters with a greater variety of core properties, focuses on 
the entropy profiles of classic cooling-flow clusters. In §2 we present our data analysis 
methods and techniques, and §3 presents our results. We flnd that the entropy proflles 
for the nine cooling-flow clusters in our sample arc remarkably similar, suggesting that the 
process preventing condensation somehow maintains the entropy profiles in a quasi-steady 
state. In §4 we discuss those results in the context of some recent theoretical models, and §5 
summarizes the paper. Throughout the analysis, we assume a cosmology with Hq = 70 km 
s^^ Mpc~^, fljn = 0.3, and fl\ = 0.7. The derived data products for the Chandra Cluster 
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Entropy Library such as the X-ray spectra, associated response files, and surface brightness 
profiles will be available through the NASA High Energy Space Archive (HEASARC), in the 
Chandra section of W3Browse^. 



2. Data Analysis 

The goal of our analysis was to derive entropy profiles as a function of radius for the 
nine clusters in our sample, with entropy quantified in terms of the adiabatic constant 

K = kTue ■ Thus, we needed to determine the electron density Ue and the gas temperature 
T as functions of radius. To do that, we first processed and cleaned the data as described 
in §2.1. Then, we divided each cluster into concentric annuli and extracted a spectrum 
in each annulus as described in §2.2. We describe our analysis of the projected spectra 
from gas within these annuli in §2.3. In order to obtain gas temperature as a function 
of physical rather than projected radius, we then performed a deprojection analysis of the 
projected spectra as described in §2.4. The deprojected cluster temperatures are not much 
different from the projected gas temperature for these clusters, presumably because the 
radial emission profiles are so steep. The resulting temperature gradients are rather coarse, 
with the number of radial bins ranging from four to thirteen. Because we desired a more 
finely-grained representation of the density gradients, we determined them by deprojecting 
the exposure-corrected source counts in the 0.5-2.0 kcV band as described in §2.5. In all the 
analysis, we use only the observations by the Advanced CCD Imaging Spectrometer (ACIS) 
backside illuminated S3 chip, which generally constrained the analysis to the inner ~ 100 kpc 
for these nearby clusters. 

2.1. Data Processing 

We reprocessed the Level 1 events files from the Chandra archive with CIAO 3.2.2 and 
CALDB 3.1.0 to obtain processed (Level 2) events files. We applied updated gain maps and 
standard grade filtering to these events. Because emission from low-redshift clusters usually 
fills the whole S3 detector, we constructed matching blank-sky background event files for 
each observation, using Maxim Markevitch's blank sky background database. We cleaned 
the cluster data of flare contamination, following Maxim Markevitch's cookbook, in order 
to match his background maps. This process involved inspecting a light curve for events 
between 2.5-7.0 keV and removing time intervals with significant background fiares, i.e.. 
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peaks with count rates > 3 cr and/or a factor of > 1.2 off the mean background level of 
the observations, using Maxim Markevitch's light curve cleaning routine lc_clean. sl.^ We 
then created a final events file for the purposes of studying the extended cluster emission by 
excluding any bright point sources near the cluster. The point sources were identified by a 
visual inspection of the image. 

2.2. Spectral Extraction 

In each cluster, wc selected a region for spectral extraction by centroiding the cluster 
emission and placing it at the center of a circular aperture with the largest radius that 
would fit entirely within the chip. Because many clusters are centered at the aim point of 
the S3 chip rather than the center of the chip, that maximum radius can be as small as ~ 2' 
(as opposed to ~ 4' for the center of the chip). For many of the clusters, this is < O.lrsoo, 
where rsoo (the radius at which the average cluster density is 500 times the critical density) is 
roughly the virial radius of the cluster. We report the maximum radius in units of arcminutes 
and kiloparsecs in Table 1. 

We then divided this aperture into concentric annuli, determining the annular bound- 
aries from background-subtracted, cumulative count profiles (0.3-8.0 keV) by fixing the num- 
ber of counts in each annulus to be at least 10,000 - 20,000 counts. The minimum number 
of total counts was chosen to allow us to test whether and where two-temperature spectral 
fits might be preferred over single temperature fits, and to allow for a simple two-parameter 
fit to the temperature gradient. 

We then extracted spectra for each annulus, along with corresponding blank-sky back- 
ground spectra, using the acisspec script. Two detector response files must be created for 
each individual spectrum, an Redistribution Matrix File (RMF) and an Ancillary Response 
File (ARF), which take into account the variability and the spatial dependencies across the 
detector. As of CIAO 3.2 and CALDB 3.0, the CIAO tool mkacisrmf must be run after spec- 
tra are extracted for each cluster using the acisspec script. The purpose of mkacisrmf is to 
generate an RMF using a CCD spectral response with two components. One component 
does not include charge-transfer inefficiency (CTI) effects; the second incorporates the spa- 
tial variation in the chip response caused by CTI. The weight map created by acisspec is 
utilized by mkacisrmf to make a count- weighted RMF. Finally the weighted ARF is created 
with mkwarf using the same weight map and the new weighted RMF, in order to match the 
energy grid of the weighted RMF. (The last step utilizing mkwarf is only necessary because 
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we fit our data using XSPEC.) For these older data on extended sources from the ACIS S3 
chip, this new procedure results in temperatures and metallicities only marginally different 
from that obtained with CALDB 2.29, CIAO 3.0. Therefore, one can be confident that, for 
the intents of this paper, the Chandra calibration has stabilized for archival data and our 
results are robust to the details of how the RMF or ARE were computed. 

However, in our experience since the beginning of this project, the evolution of the 
Chandra calibration software and the associated calibration data files has occasionally led to 
significant changes in the best-fitting spectral parameters for these same observations since 
their original publication. Wc discuss those changes in § 3.3. As of late 2005, such changes 
appear to be a thing of the past, as the cahbration of the ACIS S3 detector seems to be 
converging to consistent results. 

For the spectral fitting, wc grouped the counts into bins with a minimum of 25 counts 
per bin, and restricted the fit to 0.7-7.0 keV. 

2.3. Projected X-ray Spectra 

Once we had extracted the spectra, we fit them with plasma emission models using 
the software package XSPEC 11.3.1 for several classes of input assumptions, including one- 
temperature (1-T) and two-temperature (2-T) MEKAL models. We fit each annular spec- 
trum between 0.7 — 7.0 keV individually, so that all free parameters are independent in each 
annulus. Earlier versions of the Chandra calibration seemed to suggest that two-temperature 
models were required. However, after the Chandra calibration improved, we discovered that 
two-temperature models were no longer required to fit the vast majority of the cluster spec- 
tra in our sample. We also explored the effect of fixing the hydrogen column density (A^h) 
of X-ray absorbing gas to the Galactic value in the direction of the cluster, finding that 
discrepancies between the best-fit hydrogen column density and the assumed Galactic value 
resulted in a systematic shift in the derived cluster temperature. Again, with earlier versions 
of the calibration data and software, this source of uncertainty affects the hotter clusters 
more than the moderate temperature clusters. As the calibration improved over the years 
we worked on this project, the impact of hydrogen column density anomalies on the best- fit 
temperature diminished. As of the last two calibration versions, the best-fit A^h is consis- 
tent with Galactic A^h from Dickey & Lockman (1990). Therefore, our spectral results are 
reported for fits for which A^h was constrained to an assumed Galactic value. 

Each one-temperature fit takes four parameters: a temperature T, a heavy-element 
abundance fraction Z relative to the solar values, an absorbing column density Ah, and a 
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spectral normalization defined to be 



10-14 



/ 



HenndV 



(1) 



where rif. and nn are the number densities of electrons and hydrogen nuclei, respectively, 
in units of cm^^, and dp^ is the angular-size distance in cm. Table 2 gives the best fits for 
each annulus when A'h is fixed to the Galactic value toward the cluster taken from Dickey 
& Lockman (1990). The columns in Table 2 are as follows: 

1. Cluster name 

2. Outer radius of annulus in arcminutes 

3. A^H in units of 10^° cm'^ 

4. Best-fitting temperature in keV with 90% confidence limits 

5. Heavy-element abundance Z in solar units along with 90% confidence limits [The solar 
value of Fe/H is taken to be 4.68 x 10"^ (Anders & Grevesse 1989).] 

6. Spectral normalization N of the MEKAL model fit 

7. Statistical uncertainty of the spectral normalization N with 90% confidence limits 

8. Reduced ^ value of the fit 

9. The number of degrees of freedom in the fit 

In some cases, the values of the 1-T fits to the projected data are poor. A poor 
fit could mean either that multiple temperature components are present at a given radius 
or that the projected spectrum in that annulus comes from gas of different temperatures 
because of a radial temperature gradient. To resolve this issue, we deprojected the annular 
spectra for each cluster. 



In a cluster with a centrally peaked surface-brightness profile, the spectrum from an 
annulus of a given projected radius is dominated by gas at radii similar to the projected 
radius but includes cluster emission along the entire line of sight through the cluster. In 
order to recover the properties of the gas as a function of physical radius, one must correct 
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for this projected emission. To deproject spectra, one starts at the outermost projected 
annulus, fits an emission model and then iteratively removes the contribution of the outer 
layers from inner annuli. This procedure requires, at minimum, some sort of assumption 
about the symmetry of the cluster. For this work, we assumed that the cluster emission is 
spherically symmetric. 

To deproject the cluster data, wc fit the spectra using the built-in dcprojcction model 
projct in the software package XSPEC, tested and validated by Johnstone ct al. (2005). 
The fit simultaneously includes the projected spectra from all annuli in the cluster, and the 
initial guesses to the fit parameters were taken from the fits to the projected spectra. Table 3 
lists our best-fit parameters when iVn is fixed at the Galactic value. The columns in this 
table are the same as for Tables 2, except that the normalization in column (6) is that for 
a spherical shell, not for the individual spectrum at that annulus. The values quoted 
represent a single fit over all the spectra for the cluster instead of those for the individual 
annuli. 

The deprojected fits occasionally showed signs of instability in which the temperature 
and metallicity would oscillate, with large uncertainties, from shell to shell. Since an ac- 
curate measurement of metallicity requires more counts than an accurate measurement of 
temperature, we tied the metallicities of neighboring annuli together. The metallicity was 

constrained to be equal across groups of two or three annuli to make a common metallicity 
estimate at lower spatial resolution than temperature or normalization. Even so, the best fit 
projected temperatures occasionally exhibited instability from shell to shell. However, the 
excursions were smaller than the statistical uncertainty in the temperatures. 

We also noticed that deprojections of spectra with background corrections based on 
the deep field resulted in the outermost bin having somewhat higher normalizations (see 
also Johnstone et al. (2005)). The outermost bin in most cases is contaminated by emission 
originating even farther from the center of the cluster. To estimate the effect of cluster 
emission outside the outermost annuli, we deprojected the spectral datasets for each cluster 
where the spectrum from the outermost annulus was background-corrected based on local 
background instead of the deep fields. This choice reduced the best-fit normahzation in that 
bin but did not change the best-fit temperature in that bin. Our results in this paper are 
therefore not sensitive to the treatment of this outermost bin. We report the deprojection 
results for deep-background subtracted data only in Table 3, except for the case of Abell 1795 
which had a higher local background than average. (The temperature for the second most 
outer annulus had no upper limit for A1795 data for which all annuli had a deep-background 
correction.) The outer temperature was statistically the same for either method. 

Using this deprojection algorithm, we obtained reasonable values for fits with a 
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single-temperature plasma within each spherical shell. We tested the effect of including a 
second temperature component, and occasionally a second component in the central sphere 
improved the somewhat, but not significantly. The three clusters for which there may be 
evidence in our data for a second component are 2A0335+096, Abell 252, and Abell 2052. We 
report the results of a dcprojection analysis where we allowed a second thermal component 
in the innermost bin in Table 4. The metallicity of this inner sphere was constrained to be 
the same in both components. 

Wc did not require second components in as many cases as previous workers, particularly 
those who included soft energy bins in their fit {E < 0.7 — 0.8 kcV). This discrepancy is 
due to the immaturity of the Chandra calibration early in the mission. We found over the 
history of our own analysis that the need for a second temperature component decreased as 
the CALDB version number increased. We will discuss the sensitivity of our derived entropy 
profiles to the presence or absence of a soft component in the central cores in §3.2. 

2.5. Electron Densities from Deprojected Surface Brightness Profiles 

Electron density profiles can be derived with much higher resolution than the tempera- 
ture and abundance profiles given in Tables 2-3 because the count rate in a limited bandpass 
is much more sensitive to electron density that it is to temperature. We therefore derived 
high-resolution density profiles for each cluster by dividing our apertures into annuli of 2.5 
and 5.0 arcseconds, and using interpolated normalization-to-count-rate ratios to solve for 
the electron density corresponding to a given projected count rate. 

We created differential surface-brightness profiles from 0.5-2.0 keV to minimize their 
dependence on temperature. We corrected the counts in each annular bin for vignetting and 
small variations in the net exposure time for each bin by extracting a exposure profile from 
a normalized exposure map created assuming a mono-energetic photon spectrum of 1 keV. 
These corrections were typically less than 5% per bin. Comparisons with maps made for 0.5 
and 1.5 keV photons showed that, the systematic uncertainty induced by assuming mono- 
energetic photons is negligible. We found that the ratio between the spectral normalization 
quantity N and the count rate in this bandpass (0.5-2.0 keV) was relatively insensitive to 
temperature, changing by about 10-15% over the full range of temperature in a given cluster. 
The conversion was actually more sensitive to metallicity than temperature in this bandpass. 
Therefore, the derived electron density profiles have only limited sensitivity to the details of 
how we interpolate cluster temperatures or even how the cluster temperatures themselves 
were determined. 
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A deprojected emission profile (count rate per unit volume) was then computed for each 
cluster using a standard technique (Kriss et al. 1983). The deprojected emission profile was 
converted into electron densities using the appropriate ratio of a MekaL spectral normal- 
ization to spectral count rate in the same bandpass as the surface brightness profile. This 
conversion takes into account both temperature and abundance variations that somewhat 
affect the emissivity of the gas in the 0.5-2.0 keV energy range. 



Here we describe the results of our data analysis. We find that the clusters in our 
sample have rising temperature gradients and declining abundance gradients, in agreement 
with previous work. Wc briefly discuss those results in §3.1. Then, in §3.2, wc present high- 
resolution entropy profiles for our clusters derived by interpolating the observed temperature 
gradients onto our density-profile bins. These entropy profiles turn out to show a striking 
regularity that is best fit with a power law in radius plus a constant entropy pedestal. 
We conclude the section with a cluster-by-cluster comparison of our results with previous 
observations of these clusters. 



All the clusters in our sample have temperature profiles that rise with radius. In order to 
quantify those temperature gradients, we fit the deprojected temperature data for the seven 
clusters that had at least four annuli with a minimum of 20,000 counts with the power-law 
relation 



where Tioo is the temperature in keV at 100 h^Q kpc and rmid is the midpoint between the 
inner and outer radii of the annulus. We report the results of these fits in Table 5. This 
simple form gave adequate fits for all clusters except Abell 2052. We find that the typical 
power-law index for these temperature profiles is cct ~ 1/3 inside 4', in agreement with the 
results of Voigt & Fabian (2004). 

The abundance measurements in these clusters, which are dominated by the iron lines in 
the spectrum, generally decline with radius. The innermost regions tend to have Fe/H equal 
to 50-100% the solar value, decreasing to 30% solar at > 100 kpc. This finding generally 
agrees with previously observed trends in coohng-fiow clusters (e.g., De Grandi & Molendi 
2001). In §3.3 we discuss our abundance results for individual clusters in more detail. 



3. 



Results 



3.1. Abundance and Temperature Profiles 




(2) 
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3.2. Entropy Profiles 



Our primary new results concern the entropy gradients of cooling-flow clusters. Using 
the spectral fitting and deprojection results of §2, we derived a radial entropy profile for each 
cluster under the assumption that the temperature and density distributions were spherically 
symmetric. Strict spherical symmetry is obviously an idealization of reality that applies 
better to some clusters than to others, but we wanted to derive entropy profiles for all 
the clusters in a uniform way. To obtain the entropy profiles, we linearly interpolated 
temperature profiles onto the fine radial grid used for the surface-brightness deprojection 
in §2.5. That gave us ne(r) and T(r), from which we constructed K{r) = kT{r)ne'^^^{r). In 
order to test the effects of different binning schemes, we used bin widths of both 2.5" and 5" 
per annulus. Since deprojecting the spectra often did not give significantly different results 
from simpler analysis of projected spectra, and since projected temperature profiles were 
better behaved, we used the projected temperatures for the profile fits reported in Tables 6 
and 7. To show how much of an effect this choice had on our results, we also include the 
results for entropy profiles obtained by using the power-law fits to T(r) from deprojected 
spectra (Table 8). 

Because interpolation of temperature within the innermost temperature bin is not well 
constrained, we treated the temperature gradient in a few ways, in order to probe the 
sensitivity of our results to how we handle the temperature profile. Method 1 modeled the 
temperature gradient in the innermost bin with a linear extrapolation from the adjacent 
bin. Method 2 used a constant temperature within the innermost bin. This assumption 
may have the effect of inducing a core in the entropy profile. The differences between the 
two methods provide an estimate of the systematic uncertainties in the modelling of the 
temperature structure in the innermost bin. Method 3, where we used the best-fit power 
law to the deprojected temperatures, provides yet another estimate which also minimizes 
the central entropy core since the temperature at r = is forced to be T = in such fits. 
This third model should be thought of as a providing a lower limit on the central entropy, as 
there is no spectroscopic evidence for significant amounts of gas with such low temperatures 
in the centers of clusters. 

We fit the entropy profiles derived using Methods 1-3 with two different formulae. The 
first formula assumes that the entropy profile is a power law in radius plus a constant entropy 
pedestal Kq: 




(3) 
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The second formula sets Kq — 0, making the entropy profile a pure power law: 

Table 6 reports the best fits obtained with these two formulae using surface-brightness bins 
5" in width. Table 7 reports the best fits using surface-brightness bins 2.5" in width. Table 8 
reports the best fits using surface- brightness bins 2.5" in width and the best fit to the 
deprojected temperature profile. Our tables also report the range (in Mpc) of the fitted 
region. 

We used 1000 Monte-Carlo bootstrap reproductions of the original surface brightness 
data to quantify the statistical uncertainties of the deprojection process. The outermost de- 
projected bins in our fit are noisy and were excluded from the fit. For one cluster, Abell 2052, 
the central surface brightness deprojection is strongly influenced by the presence of cavities 
and edges. The overall fits for this cluster consequently have large values because of 
these non-axisymmetric structures, which make the azimuthally averaged surface-brightness 
profile and therefore the deprojected density profile non-monotonic at r ~ 10 kpc. 

Our entropy-profile fitting revealed a remarkable uniformity among the clusters in our 
sample, which can be seen in Figures 1 and 2, despite the variety of structures that can be seen 
in the cluster morphology. All of them have approximately the same entropy normalization 
at ~ 100 kpc: Kiqq ^ 150 ± 40 kcV cm^, Kiqq ^ 150 ± 50 keV cm^, and f^ioo ^ 140 ± 30 
keV cm^ for methods 1, 2, and 3 respectively. The best fits are generally obtained with a 
non-zero value for the entropy pedestal: ~ 7 ± 4 keV cm^, Kq ?a 11 ± 5 keV cm^, and 
~ 6 ± 4 keV cm^ for methods 1,2, and 3 respectively. They also have similar power- law 
slopes in the 10-100 kpc range: a = 1.2 ± 0.2 for fits with Kq ^ 0, and a = 1.0 ± 0.2 
for fits with — 0. Only Abell 2029 has an entropy profile marginally consistent with a 
vanishing entropy value at r = 0, having = 3.0 ±2.1 (for 5" bins) when Method 1 is 
used for temperature interpolation. These results are insensitive to the binning procedure, 
as one can see by comparing Table 6 with Table 7. The main effect of using a deprojected 
temperature profile instead of a projected temperature profile is generally (but not always) 
a slightly lower central entropy quantity Kq (compare Table 7 with Table 8.) The largest 
uncertainties come from the statistical uncertainty of the gas temperature and from the 
treatment of the temperature profiles in the central bin. 

Our results on central entropy values [Kq) pertain to the component that fills the 
majority of the volume in the central bin. In many cases, there is clearly gas of very low 
entropy {K ~ 10~^ keVcm^ near the center of the cluster in the form of Ha-emitting nebular 
filaments. Another example of cool, multiphase gas can be seen in the central 10 kpc of M87. 
The X-ray filaments along the radio source consist of gas at ~ 1 keV surrounded by a 2 keV 
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cluster atmosphere (e.g., Sparks et aL 2004), but those filaments do not constitute a large 
fraction of the volume in the 10 kpc sphere surrounding M87. 

In order to evaluate how the presence of a cool component affects the entropy derived 
for the hot volume- filling component, one can consider the emission from a two-component 
plasma with temperatures Th and Tc and normalizations Nh and Nc for the hot and cool 
components, respectively. If the two components are in pressure equilibrium, then the elec- 
tron density of the hot component is proptoN^J'^ * [1 + {Nc/Nh){Tc/Th)'^]'^^'^ . In the three 
clusters for which there is some evidence for a second temperature component, that compo- 
nent is present only within the central bin and has Nc ~ 0.5A^/,, and T,. ~ 0.5Th. Such a cool 
component comprises only ~ 10% of the X-ray emitting gas mass, and failing to account 
for it leads to a derived value of electron density that is only a few percent different from 
the actual value, if Nh is properly measured. A larger uncertainty arises in the temperature 
measurement, because the best-fitting temperature in a single-temperature model will be 
a weighted mean of Tc and T^, that depends on the spectral band of the fit. Comparing 
our single-temperature fits with our two-temperature fits shows that Th is generally ~ 50% 
higher than the best-fitting single temperature. In addition, if a cool component is present, 
then the normalizations derived from single-temperature models are overestimates of the 
true normalization of the hot component in the central bin. In the most extreme case, in 
which ~ 50 of the counts come from the cool component, the electron density inferred from 
single-temperature models would be ~ 40% larger than the actual value in the hot phase. 
Combining all these effects would mean that the central entropy of the hot component could 
be as much as ~ 80% larger than what we estimate from our single temperature models, 
if a second component of cooler gas is indeed affecting our temperature measurements, and 
the entropy of a cool component with Tc ~ 0.5Th would be ~ 50% of the single-temperature 
entropy estimate. 

Another source of uncertainty in entropy values derived under the assumption of a single- 
component plasma crops up in cases where radio plasma creates X-ray cavities. In such cases, 
the volume of the X-ray emitting plasma is smaller than assumed, meaning that the electron 
density is underestimated and the entropy is overestimated. However, this effect is not large 
compared with other sources of uncertainty. If the radio plasma displaces as much as 25% 
of the gas in a given annulus, then the true electron density would be (0.75)~^/^ ~ 1.15 of 
what was estimated, and the actual entropy would be (0.75)^/^ ~ 0.91 times the estimated 
value. 

In summary, the entropy profiles we derive here are for the X-ray emitting component 
that fills the majority of the volume at each radius. Because the temperature structure of the 
central bin is difficult to establish, the systematic uncertainties in our central entropy mea- 
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surements are approximately a factor of two. In three of our clusters, spectral deprojection 
suggests that there may be a second temperature component in the central bin at roughly 
half the temperature of the volume-filling component. Single-temperature modeling of the 
central bin could therefore be underestimating the entropy of the volume- filling component. 
Our finding that the entropy of the volume-filling component approaches a minimum value 
~ lOkeVcm^ at small radii is therefore robust to the presence of a cooler component. The 
mass of gas in a cooler X-ray emitting component must be considerably less than that in 
the hotter component. Otherwise, it would emit very bright soft X-ray emission lines that 
would show up in our spectra as a pronounced soft excess. High-resolution spectroscopy of 
the cores of many of these clusters also severely limits the amount of gas cooler than about 
1/2 — 1/3 the virial temperature (Peterson et al. 2003). In the future, it will be fruitful to 
compare these results with a parallel analysis of the predicted X-ray emission of simulated 
3D clusters with realistic temperature profiles and emission structures. We are pursuing that 
work separately. 



3.3. Individual Clusters 

The following paragraphs discuss the individual clusters in our sample. In each case, 
we give the radio and optical emission line characteristics of the cluster, which can affect 
the X-ray morphology. Then we discuss how our results on the individual clusters compare 
with previously published Chandra analyses. We also compare our results with published 
XMM-Newton results. All clusters in the sample have a central radio source in a single 
central bright galaxy; all but one (Abell 2029) has an optical emission line nebula in the 
brightest central galaxy. We also note that all but one of the clusters in our sample with Ha 
emission also have had confirmed detections of vibrationally-excited molecular hydrogen at 
2 microns (Edge et al. 2002; Donahue et al. 2000; Jaffe et al. 2001; Falcke et al. 1998; Jaffc 
& Bremer 1997), except for 2A0335+096, which has a CO detection by Edge (2001), and 
except for Abell 133, which hasn't been observed for these lines recently. 

For the most part, our results on the temperature and metallicity profiles of these clus- 
ters seem relatively robust to calibration and telescope differences. The main discrepancies 
we find stem from updates to the Chandra calibration that have reduced the prominence 
of second temperature components or anomalous absorbing column densities suggested or 
required by the data in previous analyses. For a similar reason, single temperatures derived 
from fits that included all or part of the 0.3-0.7 keV bandpass in those analyses of Chandra 
data were occasionally different from what we are now obtaining. 

Our results are in rough agreement with recently published XMM-Newton temperature 
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profiles for individual clusters. XMM-Newton observations of six of the clusters in our sample 

were analyzed by Kaastra et al. (2004) and Tamura et al. (2004): 2A 0335+096, Abcll 262, 
Abell 2052, Hydra A, Abell 496, and Abell 1795 for the first paper and 2 A 0335+096, Abell 
262, and Abell 2052 for the second. For single-temperature projected fits, we find that 
the Chandra temperatures are consistently about 0.5 keV higher than XMM temperatures 
outside 0.5'. For Abell 1795, the discrepancy is larger, ~ +1 keV, over the entire region. 
For single temperature deprojected fits, we have no discrepancy with Tamura et al. (2004) 
(Abell 262, Abell 2052, and 2A0335+096), except for the radii of 0.5 - V for 2A0335+096, 
where we obtain a temperature about 0.4 keV higher. For the 2T deprojected fits where a 
cool component is posited for the center of the cluster, we agree with Tamura et al. (2004), 
but for a ~ +0.2 keV discrepancy in Abell 262's 0.5 — 1.0' annulus. We shall show in 
the following discussion of individual clusters that our analyses rarely disagrees with that 
of other Chandra observers, so this discrepancy is part of a systematic difference between 
Chandra and XMM that appears to have reduced over time. The metallicities we find agree 
with those found by Tamura et al. (2004) for all three clusters in their sample. 

The slope and normalization of our entropy measurements are consistent with recent 
measurements from XMM. Piffaretti et al. (2005) derive and plot entropy profiles from XMM 
data for four clusters in our sample: 2A0335+096, A262, A2052, and A1795. After converting 
for different Hubble constants, and assuming that their outermost point is coincident with 
Rout listed in their Table 1, we find excellent agreement between the estimated entropy values 
plotted in their Figure 5 with our entropy measurements for those clusters. Encouragingly, 
their entropy profiles at large radii interpolate in to the profiles we measure with very similar 
slope, and the entropy measurements at overlapping radii agree in their normalization. Our 
Chandra profiles have higher resolution and but span a smaller range of radii. (Abcll 496 
is also in their sample, but we could not unambiguously identify its points in their entropy 
plot.) 

The fact that our Chandra results generally agree with the Tamura et al. (2004) analysis, 
and the entropy profiles of Piffaretti et al. (2005), and yet have a nearly constant offset in 
temperature with respect to Kaastra et al. (2004), suggests that the XMM results do not 
agree with each other. The discrepancy may arise from using different energy ranges in the 
fitting (Kaastra et al. (2004) fit 0.2-10.0 keV for every spectrum, while Tamura et al. (2004) 
fit 0.3-8.0 keV for the MOS detectors and 0.7-7.0 keV for the pn detector, similar to the 
energy range we use to fit Chandra data (0.7-7.0 keV). Furthermore, Tamura et al. (2004) 
use a more recent version of the XMM cahbration and software (SAS 5.3.0 vs SAS 5.3.3). 
The main difference between these versions is a new tool (especget) for the extraction of 
spectra and the computation of responses (RMFs and ARFs). Our consistency with later 
XMM and Chandra results suggests that the the calibrations of the two missions are leading 
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to increasingly consistent results, which is encouraging news. 

We discuss specific comparisons to pubhshed analyses of Chandra observations below. 

3.3.1. 2A 0335+096 

The cluster 2A0335+096 is poor but popular. Sarazin et al. (1995) presented VLA 
images showing that this cluster contains a double-lobed source with a nucleus. Each lobe lies 
about 12" on either side of the nucleus of the central galaxy. That galaxy has a spectacular, 
filamentary emission-line system with a bar and filaments in the central 20 kpc (Romanishin 
& Hintzen 1988). The X-ray images show two prominent cavities along with filamentary 
structures that are not correlated with the radio or optical line emission (Mazzotta et al. 
2003). The structure of the radio source suggests that it was produced by multiple outbursts. 

Chandra data on this cluster has previously been analyzed by two groups, Kawano et al. 
(2003) and Mazzotta et al. (2003). Their results disagree somewhat with each other. Our 
results agree more with the latter, who may have used a later CALDB version than the 
former group. 

Our comparison with the Kawano et al. (2003) analysis shows that our projected 1-T fit 
agrees in the region < r < 20 kpc, with temperatures between 1.65 — 1.76 keV. However, 
the temperature gradient we derive rises more quickly over the 20 < r < 120 kpc range. 
At r ^ 40 kpc we find a best-fit temperature value T = 3.0 ± 0.1 kcV while Kawano ct al. 
(2003) find T = 1.9±^om keV. At r 54 - 70 kpc our best-fit temperature is 3.8 ± 0.15 keV 
while Kawano et al. (2003) obtain 2.5 ±0.1 kcV. Our overall temperature fits arc also more 
consistent with the ASCA temperature of 2.86 ± 0.02 kcV (Horner 2001). The difi^crcncc 
between our best-fit temperatures and those of Kawano et al. (2003) arises from the version 
of the CXC Cahbration Database (CALDB) used to calibrate the data and the energy range 
used for spectral fitting. We used CALDB 3.0 in our work, while Kawano et al. (2003) used 
CALDB 2.9. Also, we fit the spectral range 0.7 - 7.0 keV, whereas Kawano et al. (2003) fit 
the 0.5 — 10 keV range. In our experience, the earliest results from Chandra for the coolest 
clusters are the most affected by improvements in the calibration. Also, including bins at 
the highest energies without many counts in them can lead to unreliable results. 

We also compared our projected and deprojected one-temperature MEKAL fit results 
with those of Mazzotta ct al. (2003) who used the Kawano et al. (2003) data. They do not 
report the version of CALDB that they used. The temperature gradient we measure has a 
similar slope to the one they obtained. Over the region of 0-240", Mazzotta et al. (2003) 
their best-fit temperature rises from 1.8 keV to 4.2 keV with a fiattening at 100" < R < 200". 
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The metallicity and A^h profiles are also similar in value and behavior. A second temperature 
component in the central region marginally improved our deprojected fit (Table 4), but the 
lower limits on the normalization of this component and the modest improvement of are 
not indicative of a secure detection, particularly since the uncertainties of the calibration are 
probably highest at the energies of interest here. 

3.3.2. Abell 133 

The center of Abell 133 hosts an impressive radio relic source spanning 55 kpc. Slee 
et al. (2001) conclude from 4" resolution VLA observations that the central galaxy is not 
the current source of the relic but may have been where the radio source began. ROSAT 
X-ray observations show that the radio source is clearly interacting with the ICM (Rizza 
et al. 2000). The central galaxy also has a compact, low-ionization emission-line source (Hu 
et al. 1985). 

Fujita et al. (2002) report Chandra observations of Abell 133. We find statistical agree- 
ment with their projected temperature, metallicity, and A^h profiles for the fits in which 
Fujita et al. (2002) restrict the spectral fitting to photons with energies greater than 0.9 
keV. However, their fits for data which include the less energetic photons disagree with ours. 
This finding is in accord with the general trend that early Chandra spectral fitting results 
involving soft X-ray data are probably not reliable because the low-energy calibration for 
Chandra was not well characterized until at least 2004. Hence, the energy range chosen for 
any given fit can have a significant impact on the best-fit temperatures, absorption column 
A^H) arid metallicity values. The fact that the fits of Fujita et al. (2002) that were restricted 
to higher energies agree with ours supports our suspicion that the Chandra calibration of 
the soft X-ray bandpass has improved with time. 

3.3.3. Abell 262 

Abell 262 hosts the weakest radio source of our sample, the doubled-lobed source B2 
0149+35 (Parma et al. 1986; Fanti et al. 1986), but it also has an impressive emission line 
nebula (Plana et al. 1998). Blanton et al. (2004) find that the radio source is anti-correlatcd 
with the X-ray emission in their Chandra observations and but it is correlated with optical 
([N II]) emission. 

In order to have enough bins to fit a 2-parameter power law to the deprojected tem- 
perature profile, the spectral fits reported for A262 have 10,000 counts per spectrum instead 
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of 20,000. Our derived temperatures are consistent with the analyses of the Chandra data 
in Blanton et al. (2004). Our projected one-temperature fits with A^h fixed to the Galactic 
value did not describe the inner annulus between and 0.35' (< 7 kpc) well, with a reduced 
of 2.44. But, as Blanton et al. (2004) also found, adding a second temperature component 
to the projected data improved our fit with a reduced of 1-06 for 107 degrees of freedom. 
We also find good general agreement between their mean temperatures and metallicity and 
the average temperature and metallicity of our radial fits. 

Fits to deprojected spectra arc typically better than those for inner projected spectra, 
suggesting that departures from singlc-tcmpcraturc emission in the projected spectra stem 
primarily from a steep radial temperature gradient, which leads to superpositions of multi- 
temperature plasma along the line of sight. This trend was also true for Abell 262. In the 
central bin for Abell 262 (i.e., < 0.35'), adding a second temperature component to the 
deprojected model only marginally improved the fits to the projected spectra (Table 4). In 
these deprojected spectral fits, a single-temperature model in each radial shell gave a reduced 
■y^ — 1.28 for 923 degrees of freedom, while a two-temperature plasma model in the inner 
two shells give — 1.22 for 921 degrees of freedom. 

3.34. Abell 496 

A compact radio source (smaller than ~ 1.5" based on the beam size reported) inhabits 
the central galaxy in Abell 496 (O'Dea et al. 1995). The central galaxy is also the locale for 
a bright emission-line nebula (e.g. Heckman et al. 1989). Since the observation of Abell 
496 had only 60,000 total counts, we first divided this cluster into three regions. When 
our spectral analysis of these three regions turned out to be well-behaved, we expanded 
the analysis to six regions of 10,000 counts each in order to fit a 2-parameter power- law 
temperature profile, and it is this result which is reported in our spectral fits. 

Dupke & White (2003) present the original Chandra analysis of these data. Comparison 
of our projected one-temperature fits with those from their analysis show good agreement 
across the region — 150 kpc. Our temperature and metallicity profiles are also consistent 
with the results of deprojected, 1-T fits to the XMM-Newton data (Tamura et al. 2001a). 

3.3.5. Abell 1795 

The central galaxy of Abell 1795 contains a famous extended filament of optical line 
emission, mapped using long-slit spectroscopy by Cowie et al. (1983). The central bright 
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galaxy also hosts a compact radio source 4C 26.42 (Caswell Sz Wills 1967). Chandra obser- 
vations by Fabian et al. (2001) revealed a 40" X-ray filament that substantially overlaps the 
optical emission- line filament. 

The results we report on this cluster agree well with previous X-ray results. Our pro- 
jected and deprojected one-temperature and two-temperature fits are in excellent agreement 
with the Chandra data analysis by Ettori et al. (2002). They also agree with the projected 
one-temperature, A^Hi and metallicity fits from XMM-Newton analysis of Abell 1795 by 
Tamura et al. (2001b). 

3.3.6. Abell 2029 

Abell 2029 occupies a unique niche in our sample because it alone has no trace of an 
optical emission line nebula, but it is still a luminous radio source. Its X-ray image is 
remarkably smooth, exhibiting little of the structure seen in some of its sister cool core 
clusters. Therefore it has been a textbook example for the quest to determine the self- 
interaction cross-section of dark matter by fitting dark matter potentials to the enclosed 
mass inferred from X-ray data (e.g., Lewis ct al. 2003). 

The work of Lewis et al. (2002, 2003) on Chandra observations of Abell 2029 adopts the 
APEC spectral model within XSPEC and fixes N// to the Galactic value. Our projected and 
deprojected temperature fits, along with our metallicities, agree with that group's analysis. 
Our temperature fits confirm the fiattening of the temperature profile inside of 18" found by 
Lewis et al. (2003). 

3.3.7. Abell 2052 

The central galaxy of Abell 2052 hosts the radio galaxy 3C 317. Venturi et al. (2004) 
found a parsec-scale bipolar radio source with a radiative age of 170 years in VLBA observa- 
tions of this galaxy and suggested that this source is a restarted radio galaxy. The kpc-scale 
appearance of the source is that of an amorphous halo with a bright core. The Chandra 
X-ray emission shows two bubbles in the ICM (Blanton et al. 2003a, 2001). This galaxy also 
is the home of a bright emission-line nebula (e.g. Heckman et al. 1989). 

Our projected temperature profiles for fixed A^h agree with the Chandra analysis of 
Blanton et al. (2003b). We also find that our best-fit metaUicities track theirs, increasing to 
0.75 solar at r «i 30, then falling to values around 0.45 solar at larger radii. Our deprojected 
temperature and metallicity profiles also agree with theirs. In particular, we also see a 
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notable increase in the deprojected metallicity at r ~ 30 kpc. 

We also obtained a minimal but likely insignificant improvement in when we included 
a second temperature component in the innermost bin for the deprojection analysis (Table 4). 
The lower limit on the normalization of this component is quite low, which suggests that the 
detection should be treated as an upper limit for the purposes of this paper. 

3.3.8. Hydra A 

Hydra A was one of the first clusters to be observed by Chandra, and it was the cluster 
that sparked our interest in doing this study. Many people have used the Hydra A data as 
comparison for their theoretical predictions; we wanted to provide a larger sample of clusters 
with similar published measurements. 

Our deprojected single-temperature fits agree with the Chandra results for Hydra A of 
McNamara et al. (2000). We also see the temperature jump at r 70 kpc with a subsequent 
decrease in temperature at r ~ 100 kpc. Likewise, comparing our results with those of David 
et al. (2001) for projected 1-T and 2-T MEKAL fits show similar agreement. Our best-fit 
values for Nj^, metallicity, and temperature are consistent with the values in that paper. 
However, we do not find an improvement in for the innermost annulus, 0" < r < 20", 
when we allow for a second spectral temperature component. This finding is consistent with 
our results on second temperature components in the cores of other clusters in our sample 
owing to the changing cahbration in the soft energy band on Chandra's ACIS-S detector. 
The entropy profile from David et al. (2001) is completely consistent in normalization and 
shape with the one we present here. 

3.3.9. PKS 0745-191 

PKS0745-191, at z = 0.1028, is the most distant cluster in our sample. It hosts a 
powerful radio source and a luminous emission hne nebula. The results of Hicks et al. 
(2002) on Chandra data for PKS 0745-191 are consistent with those we obtain from our 
projected one-temperature spectral fits. The metallicity profiles are also similar in that 
they remain nearly constant with a value of RiO.45 solar. The projected two-temperature 
fits are not relevant here because we did not find a second temperature necessary for any 
annulus. Comparison with the results of Chen et al. (2003) for XMM-Newton data for PKS 
0745-191 yields consistency between each group's results. For projected and deprojected 
one-temperature spectral fits with fixed Nh, we find a temperature profile similar to Chen 
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et al. (2003). Our metallicity profiles are also consistent in values and trend with radius. 

4. Discussion 

The main new result emerging from this Chandra study of the core entropy profiles in 
cooling-flow clusters is that the profiles are quite similar, with several interesting features 
in common. Despite the fact that these clusters range in temperature from 2.2 keV to 7.4 
keV, their core entropy profiles all have a similar normalization (fs 150 — 160keVcm^) at a 
radius of 100 kpc, they all have similar power-law slopes within that radius, and the profiles 
generally tend to flatten to a constant value ^6 — 10 keV cm^ at the centers of the clusters. 
In this section we examine these features more closely and explore their implications. 

4.1. Observations of Non-Zero Central Entropy 

The presence of a non-zero central entropy pedestal in Chandra observations of Hydra A 
was noted by David et al. (2001). However, it has not generally been recognized as a common 
feature in cooling-fiow clusters. For example, PifTaretti et al. (2005) analyzed the entropy 
profiles of thirteen cooling-fiow clusters observed with XMM-Newton, finding that they were 
adequately fit by a pure power law with a ^ 0.95 without the need for a central entropy 
pedestal. However, the lower spatial resolution of XMM-Newton (4") makes it harder to 
resolve the central ~ 10 kpc where the entropy pedestal dominates the entropy profile. Other 
Chandra studies suggest that non-zero central entropy might not be universal in cooling- flow 
clusters. One possible counterexample from our own study is Abell 2029. Another possible 
counterexample is Abell 478 (Sun et al. 2003b; Sanderson et al. 2004). Also, there are group- 
scale objects which seem to have entropy profiles that tend to zero entropy at the center 
(Sun et al. 2003a; Mahdavi et al. 2005). 

Our finding that non-zero central entropy is common, if not universal, among cooling- 
fiow clusters rests on the fiattening of the observed surface-brightness profiles within ~ 10 kpc 
of the cluster's center, which typically corresponds to an angular radius of ~ 10" in these 
low-redshift clusters (see Figure 3). This flattening implies that the electron density proflle 
docs not diverge at the center, and therefore that the central entropy tends to some minimum 
value. In order to verify that the flattening we observe is not an artifact of the deprojection 
procedure, we compared our observed surface-brightness profiles with the profile one would 
expect if the power-law behavior observed at larger radii continued to hold all the way to 
r — 0. For this comparison, we used a power-law core model with oc r'^ and T(r) oc r^/^. 
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so that K{r) oc r. Assuming a thermal bremstrahlung emissivity (oc T^/^) then yielded a 
predicted surface-brightness profile >S'(^) oc 9~^^^. (If we had included line radiation and iron 
gradients, the central surface-brightness profile would have been even more sharply peaked.) 
Wc then convolved this power-law model with a two-dimensional Gaussian to simulate the 
effects of a point-spread function. To be conservative, we simply approximated the Chandra 
ACIS-S point spread function as a simple Gaussian, e~^^^'^^ with a = 1". This assumed PSF 
is actually a little broader than the actual one, so we are slightly overestimating the effect 
of PSF smearing on these profiles. Nevertheless, the surface-brightness profile predicted by 
this simple power-law model is significantly more peaked than the actual observed surface- 
brightness profiles, even if we do not include the divergent fiux from the central pixel in 
the calculation. Figure 4 shows that comparison with both the observed surface-brightness 
profiles and the theoretical profile normalized to unity at r = 10". Therefore, unless the 
gas temperatures in the cores of these clusters drop much more rapidly than the power-law 
model we assumed (or indeed more than has been measured in high-resolution spectroscopy), 
then they typically have nearly constant central entropy values of ~ 6 — 10 keV cm^. 

One could argue that this result is not completely unexpected. For some time now, X- 
ray astronomers have been successfully fitting "beta-model" and even "double-beta model" 
profiles to X-ray surface brightness profiles of clusters (e.g., Mulchaey & Zabludoff 1998; 
Mohr et al. 1999; Xue & Wu 2000). The standard "beta-model" for a cluster's X-ray surface 
brightness I{9) is proportional to (1 -|- 9 /9c)~^^'^^^'^ , where 9c is the projected core radius 
in appropriate units (Sarazin 1988; Cavaliere & Fusco-Femiano 1976). The X-ray surface- 
brightness data for clusters and groups have consistently exhibited evidence for cores in the 
electron density profile, which is relatively insensitive to the temperature profile. The key 
feature of the beta-model is its flat central core. These high resolution Chandra observations 
show that even those clusters with significant peaks in X-ray surface brightness still have 
cores at 6* < 10". 



The observed central entropy levels in these nine clusters suggest that the heating mech- 
anism which prevents most of a cooling-flow cluster's core gas from condensing is episodic 
(see also, David et al. 2001; Kaiser & Binncy 2003). Intracluster gas of entropy K and 
temperature T that radiates pure thermal bremsstrahlung emission has a cooling time 
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The bulk of the gas currently at the centers of these X-ray clusters therefore will not begin 

to condense for at least ~ 100 Myr. That timcscalc for the introduction of feedback is 
consistent with the periodicity of AGN feedback inferred from X-ray studies of the cavities 
associated with the radio sources at the centers of cooling-flow clusters (e.g., Birzan et al. 
2004). Voit & Donahue (2005) show that an outburst of kinetic power from an AGN at 
the level of ~ 10'^^ ergs~^ can produce such an entropy pedestal through shock heating that 
raises the core gas of a cooling- flow cluster by a uniform increment of ~ 10 keV cm^. 

Many clusters have much higher central entropy levels, which correspond to cooling 
times greater than the age of the universe, meaning that they were never suspected of 
harboring cooling flows. One classic example is the Coma cluster. It is the nearest, richest 
X-ray cluster, and it has two bright central cluster galaxies rather than a single dominant 
one. The central entropy in the Coma cluster implied by the ROSAT X-reiy observations of 
Brief et al. (1992) is SAOtH^h^Q^^ keV cm . Such a large central entropy value is difficult to 
generate purely through feedback, requiring an AGN outburst > lO^^ergs"^ in the heating 
framework of Voit & Donahue (2005). We therefore suspect that Coma's core achieved its 
high entropy level through merger shocks or some other dynamical means. 

However, there are other clusters with central entropy levels intermediate between the 
Coma cluster and the cooling-flow clusters in the present sample. These came to our attention 
when we were investigating examples of cooling-flow clusters without evidence for feedback. 
Every cluster in the sample we present here has a certain amount of central AGN activity, 
indicated by the radio power from the central galaxy. Most of them (all but Abell 2029) also 
show evidence for active star formation in the form of optical emission-line nebulae. For the 
purposes of this discussion, we will call such clusters "active clusters." In order to isolate 
the relationship between AGN activity and the cooling-flow phenomenon, we used Chandra 
to observe two other clusters, Abell 1650 and Abell 2244, that had been classifled as cooling- 
flow clusters by Peres et al. (1998) but that did not contain measurable radio power from 
a central source or an optical emission-line nebula. We will call these two clusters "passive 
clusters." 

Donahue et al. (2005) show that the two passive clusters have substantially higher 
central entropy values than the active clusters. These central entropy levels, amounting to 
30-50 keV cm^, correspond to central cooling times ~ 1 Gyr, suggesting that these clusters 
show no signs of feedback because it is not currently necessary to prevent condensation and 
may not have been necessary for quite some time in the past. One possibility is that these 
clusters were heated by an extraordinarily strong AGN outburst (~ lO'^^ergs"^) that raised 
the central entropy level to > 50 keV cm^ roughly a Gyr or more ago, so that the dynamical 
traces of that feedback event have now dissipated. Another possibility is that gravitationally 
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driven effects like merger shocks (e.g., Buote et al. 2005) have kept the central entropy 
relatively high for the last several Gyr and that these clusters will become more like classic 
cooling- flow clusters about a Gyr from now (see also McGarthy et al. 2004). 

4.3. Slopes of Core Entropy Profiles 

We find a mean power-law slope for core entropy profiles in our sample of 1.2 — 1.3 ±0.2 
when the central entropy is allowed to be non-zero and 0.9 — 1.0 ± 0.2 for a pure power- law 

fit that goes to zero at the origin. Our result for a pure power law agrees with Piffaretti 
et al. (2005), who found a ~ 0.95 in XMM-Newton observations. However, the addition of a 
~ 10 keV cm^ central entropy pedestal along with a slight steepening of the power-law slope 
clearly produces a better fit to the Chandra data, as can be seen in Figure 1. 

In either case, the power-law slope is similar to the power-law index of a ~ 1.1 ob- 
served for the entropy profiles at larger radii in clusters. That index seems to be a natural 
consequence of gravitational structure formation, which naturally produces profiles having 
K{r) oc r" with a ~ 1.1 outside the cores of clusters (e.g., Pratt & Arnaud 2003). However, 
it is not clear why that power-law behavior should continue inside the core, where the physics 
of cooling and feedback ought to dominate the thermodynamics of intracluster entropy. 

4.4. Entropy Normalization at 100 kpc 

The agreement between the entropy levels at radii of 100 kpc in these nine cooling- 
fiow clusters indeed appears to be related to the non-gravitational processes that break 
self-similarity in the whole population of galaxy clusters. Gravitational structure formation, 
acting alone, should produce nearly self-similar clusters whose whose virial radii scale with 
mass as Tyir oc M^/^ and whose temperatures then scale with virial radius as T oc r^^j.. 
Because the gas density at a given fraction of the virial radius should be the same for all self- 
similar clusters, the gas entropy at that fraction of the virial radius should be K{r/r^ir) oc T. 
For a power-law slope in entropy of a = 1.1, the entropy at a given physical radius in self- 
similar clusters should then scale as K{r) oc T^-°'/'^ oc T^-^^. Yet we see no such systematic 
trend in our sample, even though the clusters range over a factor of three in temperature. 
This finding implies that non-gravitational processes play a role in regulating the entropy 
level at 100 kpc. 

That result should not be surprising, in light of the fact that cooling and feedback are 
thought to alter the luminosity-temperature relation of clusters through their effects on the 
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core entropy distribution (e.g., Evrard Sz Henry 1991; Kaiser 1991; Voit et al. 2002; Borgani 
et al. 2004). According to the model of Voit & Bryan (2001), the entropy levels at the core 
radii of clusters should be related to the entropy level at which the cooling time of the gas 
equals a Hubble time, which is 



for pure bremsstrahlung cooling. In that case, one expects entropy at the core radius of a 
cluster to be oc T^/^, and observations of entropy at a typical core radius of O.lryir do show a 
general correspondence between Kc{T) and ii'(O.lrvir) (Ponman et al. 2003; Voit & Ponman 



Such a scaling of entropy within the cores of clusters is in accord with the agreement we 
find in entropy levels at 100 kpc. If the entropy profiles within the cores of clusters really do 
scale as K{r) oc T^/^(r/rvir)^ \ then we expect ir(lOOkpc) oc T°-^^ (see also Voit & Donahue 
2005). In other words, the observed lack of a trend in temperature in the entropy level at 
100 kpc is consistent with the idea that cooling and feedback regulate the entropy levels 
within the cores of these clusters. 



The framework for cooling and feedback presented by Voit & Donahue (2005) was mo- 
tivated by these observational findings and supports the notion that episodic AGN heating 
is responsible for governing these core entropy profiles. There we show that the core entropy 
profile one expects when pure radiative cooling is unopposed by feedback forms a lower bound 
on the set of entropy profiles we observe here. That bounding profile was computed by Voit 
et al. (2002) by simply allowing radiative losses over a Hubble time to remove entropy from 
the baseline entropy profile characteristic of gravitational structure formation. The observed 
profiles converge to this bounding profile at > 100 kpc and depart from it at smaller radii, 
where the bounding profile drops to zero entropy as r — > 0. Thus, it would appear that some 
sort of feedback is preventing the observed entropy profiles from converging to the bounding 
profile. 

The magnitude of the central entropy level is an important clue to the feedback mecha- 
nism at work. We show in Voit & Donahue (2005) that adding a constant entropy pedestal 
of 10 keV cm^ to the bounding profile reproduces the behavior of the observed profiles at 
small radii. In the framework we suggest, shock heating by an AGN outflow with constant 
kinetic power naturally produces a constant entropy pedestal out to ~ 30 kpc. As mentioned 




2/3 



(6) 



2003). 



4.5. A Framework for Cooling and Feedback 
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above, a kinetic power output ~ 10 ergs" is required to raise the inner entropy level by 
~ lOkeVcm^. Because that inner entropy level corresponds to a cooling time of ~ 10^ years, 
episodic AGN outbursts on this timescale are needed to maintain the quasi-steady nature of 
the entropy profiles of these active clusters. 

5. Summary and Conclusions 

We present temperature gradients, metallicity gradients, flux normalizations, and en- 
tropy profiles for a sample of X-ray luminous, nearby clusters of galaxies. Because we selected 
this particular sample from the Chandra archive to have large numbers of X-ray counts and 
singly-peaked surface brightness profiles in order to derive high-quality core entropy pro- 
files, wc ended up with a sample of nine classic cooling-fiow clusters. All of them show 
evidence for feedback, with radio emission from the central galaxy in all cases and central 
emission-line nebulae in eight out of nine cases, with Abell 2029 being the only exception. 
We demonstrate that the entropy profiles in the cores of these clusters are quite similar, with 
power-law slopes in radius of ~ 1.0 — 1.3, fiattening to a central entropy plateau ~ 6 — 10 
keV cm^ in the central 10 kpc or so. We suggest that the core entropy levels are maintained 
by periodic feedback from a centrally located AGN, with a duty cycle of about 10^ years. We 
demonstrate that these non-zero central entropy levels are not an artifact of finite angular 
resolution of Chandra; in fact, it is the exquisite resolution of Chandra which allows us to 
unambiguously detect these central entropy levels. 

We suggest that classifying clusters of galaxies based on their central entropy levels is a 
promising way to identify the mechanisms that prevent gas from condensing at their cores. 
The "active clusters" in the sample wc present here, with entropy levels ~ lOkcVcm^, all 
show evidence for recent feedback, but that feedback has not produced entropy inversions in 
the azimuthally averaged entropy profiles. We have observed two other clusters with Chandra 
that have been classified as cooling-fiow clusters, based on the fact that their central cooling 
time is less than a Hubble time, but that show no evidence for recent feedback. Those 
"passive clusters" turn out to have central entropy levels of ~ 30 — SOkeVcm^, even though 
they are as regular in appearance as Abell 2029, which has low central entropy (Donahue et al. 
2005). Such elevated central entropy levels can be produced by an especially strong episode 
of AGN feedback sometime in the past (Voit & Donahue 2005). It is also plausible that 
these elevated central entropy levels could be preserved by thermal conduction (Donahue 
et al. 2005). Then there are objects like the Coma cluster, with central entropy levels 
~ 350keVcm^ (Briel et al. 1992), which are difficult to achieve with AGN heating and 
would therefore seem originate through gravitationally-driven processes like merger shocks. 
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A larger survey of clusters with a greater variety of central entropy levels are needed to 
explore these issues. 
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Fig. 1. — Entropy profiles, where K = kTue . The sohd (red) triangles are the profiles 
derived assuming a temperature gradient within the central temperature bin (Method 1). 
Circles show profiles derived assuming a constant temperature in the central temperature 
bin (Method 2). The solid and dotted hues are fits to the Method 1 and Method 2 profiles, 
respectively, using a power law in radius plus a constant entropy level; the short-dashed 
and long-dashed are fits to the Method 1 and Method 2 profiles, respectively, using a pure 
power-law model that falls to zero entropy at r = 0. 
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Fig. 2. — Entropy profiles, where K = kTue . The sohd (red) triangles are the entropy 
profiles where the temperature was derived by fitting a simple power law (T oc r") to the 
deprojected temperatures (Method 3). The solid lines is a fit to the Method 3 profiles using a 
power law in radius plus a constant entropy level; the long-dashed line is a fit to the Method 
3 profile using a pure power- law model that falls to zero entropy at r = 0. 
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Fig. 3. — Background-subtracted X-ray surface brightness profiles, in units of X-ray counts 
per square arcsecond, 0.5 — 2.0 keV, where the annuh are 2.5" wide. 
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Fig. 4. — The solid line is the predicted surface brightness model for bremstrahlung radiation 
from a plasma with n ~ and T ~ r^^^. The dotted lines are the surface brightness profiles 
for all 9 clusters in this paper. All surface brightness profiles have been normalized to 1 at 
r = 10". Note that the predicted model exceeds the observed profiles in the central 10". The 
cluster with the most peaked surface brightness profile is Abell 2029. The cluster with the 
odd hump at r > 10" is Abell 2052. 
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Table 1. Clusters Included in Entropy Library 
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Table 2. Results for 1-T Projected, A^h Fixed MEKAL Fits 
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1.27C-03 


8 000-05 


-^•-^-^-0.13 


1 


79 


129 




21.12 




2 03+?,-°I 


1.46C-03 


l.OOe-04 


-^•^-^ — 0.14 


1 


12 


137 




ou.^o 




9 r,c+O.OS 


1 Qfio^O'? 
1 .OOc-UO 






1 


13 






40.64 




2 49+!!'"» 


2.29e-03 


l.lOe-04 


0.86+'ojj 


1 


01 


169 




59.76 




2 00+0.06 
— (J.OG 


4.93C-03 


1.50C-04 


0.51+"-;:? 

— 0.05 


1 


38 


223 


A496 


15.57 


4.80 


9 07+0.09 
^■■J ' -0.10 


4.86C-03 


2.80C-04 


0-95^0 is 


1 


39 


138 




28.11 




3 09+!!-;f 


5.48C-03 


2.60e-04 


79+;!-lii 

— 0.12 





85 


147 




41.78 




00+0.19 
'-'■"^ — 0.19 


5.42e-03 


2.60e-04 


0.95+Q'i^ 


1 


04 


156 




57.35 




4.06 n 01 


5.93e-03 


2.50e-04 


— U. 13 


1 


23 


160 




78 24 




4 S2+0-31 


U. J ( c-uo 


z . ouc-u i 


n C9+O.I6 


1 


08 


172 




94.95 




'-'■^^-0.37 


4.92C-03 


2.40C-04 


"•'-'"-0.18 


1 


19 


155 


A1795 


20.86 


1.17 


3.56+!!t^ 


4.47C-03 


1.40C-04 


o.62+;;-"s 

— 0.09 


1 


06 


186 




33.81 




4 33+!!'}'^ 


4.61e-03 


1.30e-04 




1 


04 


197 




46.04 




4 66+!!'?i 


4.51e-03 


1.40e-04 


"• ' ^-0.11 


1 


11 


199 




59.71 






4.79C-03 


1.20C-04 


46+"-"' 

"•™-0.09 


1 


01 


206 




74.82 




5.34+0-29 

— U.ZD 


4.92C-03 


1.30C-04 


o.47i;!-i'i' 





86 


210 




93.52 




°-°'J-0.30 


4.96C-03 


1.30e-04 




1 


03 


214 




114.38 




o.oy_,) 29 


4.92e-03 


1.40e-04 




1 


19 


216 




139.56 




6 57+0'42 
"■'"-0.38 


5.32e-03 


1.30e-04 


O-31+on 


1 


02 


226 




171.22 




6-39+rj 


3.85C-03 


1.60C-04 


o.5ot;-« 


1 


11 


231 




215.82 




0-^' -0.38 


6.79C-03 


1.50C-04 


i2+"-09 

"•-'-^-0.08 


1 


09 


252 


A2029 


9.52 


3.15 


4 14+0.28 
^■-'-^-0.27 


1.28e-03 


8.00e-05 


i-32i;5i? 


1 


31 


124 




15.58 




5 04+0.58 
^■y^-0.48 


1.72e-03 


9.00e-05 


0.65lo 21 


1 


21 


138 




20.78 




6.04+;]'6? 


1.36e-03 


9.00e-05 




1 


04 


129 




34.63 




'■-'-0-0.39 


4.99C-03 


1.40C-04 


o.67i;f;i 


1 


05 


225 




48.48 




7 Q7+0.43 


5.02C-03 


1.40C-04 


o.59t;!;i^ 





97 


226 




62.34 




7 OQ + 0.51 


5.24C-03 


1.30e-04 


0.40+oii 


1 


14 


227 




77.92 




7 70+0.51 
' ■ '■^-0.46 


5.20e-03 


1.40e-04 


o-59+;;-it 





93 


226 




95.24 




c:r+0.68 
8-00_o 5g 


5.23e-03 


1.50e-04 


0.49l'^;l« 


1 


12 


229 




115.15 




« 19+0.71 
»-i^_0.61 


4.92e-03 


1.40e-04 


44+0 -15 


1 


26 


230 



-39- 



Table 2 — Continued 



Name Rout Nu Tx Norm (TNorm Z x^ed 

kpc lO^o cm-2 kcV Zq 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 





138, 


,53 




8, 


,67 


+0, 

-0, 


68 
,60 


5, 


.63c- 


■03 


1, 


.40c- 


■04 





,36 


+a. 

-0 


13 
.13 


1, 


.05 


235 




167, 


,10 




9, 


,19 


+0, 

-0, 


,79 
,66 


5, 


.71c- 


■03 


1, 


.50c- 


■04 





,38 


+ 0. 
-0 


.14 
.14 


0. 


.98 


240 




202, 


,60 




8, 


,74 


+0, 
-0, 


,71 
,60 


5, 


,85e- 


■03 


1, 


,60e- 


■04 


0, 


.44 


+0. 
-0. 


.15 
.14 


1, 


,12 


244 




259, 


,74 




9, 


,01 


+0, 

-0, 


,65 
,57 


7, 


,60e- 


■03 


1, 


,80e- 


■04 





.40 


■fO. 
-0. 


,14 
,12 


1, 


,05 


275 


A2052 


13, 


,90 


2.85 


1, 


,67 


+0, 

-0, 


,04 
,04 


2, 


,42e- 


■03 


9, 


,00e- 


■05 


0. 


.49 


■fO. 
--0 


.05 
.05 


2, 


.59 


150 




20, 


,22 




1, 


,91 


+0, 
-0, 


,05 
,05 


2, 


.54e- 


■03 


9, 


.00c- 


■05 


0. 


,59 


-fO. 
-0 


.06 
.06 


1, 


.83 


164 




26, 


,96 




2, 


,77 


+0, 
~0, 


,08 
.08 


2, 


.24c- 


■03 


8, 


.00c- 


■05 


0. 


,97 


-fO. 
-0 


.10 
.10 


1, 


.30 


181 




35, 


,38 




2, 


,93 


+0, 
-0, 


,09 
,09 


2, 


,55e- 


■03 


8, 


,00e- 


■05 


0, 


.68 


+0. 
-0. 


.08 
.07 


1, 


,17 


186 




46, 


,75 




3, 


,23 


+0, 

-0, 


,11 

,11 


2, 


,59e- 


■03 


9, 


OOe- 


■05 





.72 


+0. 
-0. 


.10 
.08 


1, 


,12 


197 




59, 


,39 




3, 


,34 


+0, 

-0 


12 

.11 


2, 


.78c- 


■03 


8, 


.00c- 


■05 


0. 


,57 


-fO. 
-0 


.08 
.08 


1, 


.08 


206 




73, 


,71 




3, 


,51 


+ 0, 
-0, 


13 

.12 


2, 


.94c- 


■03 


9, 


.OOc- 


■05 


0. 


,50 


-fO. 
-0 


.08 
.07 


1, 


.09 


209 




96, 


,88 




3, 


,33 


+0, 
-0, 


10 
.11 


4, 


.51c- 


■03 


1, 


.10c- 


■04 


0. 


,46 


-fO. 
-0 


.06 
.05 


1, 


.15 


248 


Hydra A 


19, 


,55 


4.84 


3, 


,06 


+0, 
-0, 


11 
.10 


5, 


,50e- 


■03 


1, 


,50e- 


•04 


0, 


.49 


■fO. 
-0. 


.06 
.06 


1, 


,12 


186 




33, 


,59 




3, 


,17 


+0, 

-0, 


.12 

.11 


5, 


,38e- 


■03 


1, 


,50e- 


■04 





.37 


■fO. 
-0. 


.05 
.06 


1, 


,28 


182 




49, 


,47 




3, 


,47 


+0, 

-0 


14 

.13 


5, 


.16c- 


■03 


1, 


.40c- 


■04 


0. 


,41 


+0. 
-0 


.07 
.06 


1, 


.26 


191 




72, 


,69 




3, 


,59 


+ 0, 
-0 


15 

.14 


5, 


.38c- 


■03 


1, 


.50c- 


■04 


0. 


,35 


-fO. 
-0 


.06 
.06 


1, 


.17 


195 




122, 


,16 




3, 


,66 


+0, 

~0, 


14 
.13 


8, 


.45c- 


■03 


1, 


.90c- 


■04 


0. 


,28 


-fO. 
-0 


.05 
.05 


1. 


,00 


242 


PKS0745-191 


31, 


,75 


43.39 


4, 


,08 


+0, 
-0, 


14 
.13 


1, 


,17e- 


■02 


3, 


,00e- 


■04 


0. 


.50 


+0. 
-0. 


.07 
.07 


1, 


,24 


261 




57, 


,83 




5, 


,87 


+0, 

-0, 


.28 
.26 


1, 


,12e- 


■02 


3, 


OOe- 


■04 





.40 


+0. 
-0. 


.08 
.07 


0, 


.99 


276 




91, 


,85 




6, 


,96 


+0, 

-0, 


.37 
.34 


1. 


,09c- 


■02 


2, 


.00c- 


■04 


0. 


,50 


+0. 
-0 


.09 
.09 


1. 


,03 


286 




147, 


,42 




7, 


,47 


+ 0, 

-0 


44 
3<) 


1 


,20c- 


■02 


3. 


.00c- 


■04 





,37 


-fO. 

-n 


.08 
.08 


1. 


,14 


293 




:U(), 


,21) 




s 


„')(i 


■ (1 
-0 


■Hi 
41 


2 


,()ic- 


-02 




,0()c- 


■n 1 


0, 


,2(i 




- u 


117 


i 


,0,') 


■">)."■) 
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Table 3. Results for 1-T Deprojected, A^h Fixed MEKAL Fits 



Name 
(1) 

2A0335+096 



A133 



A262 



A496 



A1795 



A2029 



Rout 


Nh 




Norm 


^ O'VTYh 


Z 




d.o.f 


kpc 


1020 cm-2 


keV 






Zq 






(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


12.44 


18.11 


1 07+0. 06 
^■•^'-0.04 


2.55e-03 


2.50C-04 


-0.08 


1.34 


1970 


19.07 


- 


i.69irj 


3.62e-03 


2.70e-04 


- 


- 


- 


25.71 


- 


i-97l°;J° 


6.00e-03 


3.50e-04 


o.87t°:S° 


- 


- 


33.17 


- 


9 90+0.13 


5.15C-03 


3.00C-04 


- 


- 


- 


42.29 


- 


^- '^-0.14 


6.36C-03 


2.60C-04 


0.86l»:»^ 


- 


- 


54.31 


- 


2 q5+0-l» 


6.51C-03 


2.50C-04 


- 


- 


- 


70.48 


- 


fjQ+o.ai 

0.00_Q 3Q 


5.77e-03 


2.30e-04 


- 


- 


- 


90.80 


- 


Q fiQ+0.23 
O.D9_Q 22 


8.58e-03 


2.30e-04 


f;t;+0.06 


- 


- 


165.84 


- 


4 4c;+0-l2 


2.50C-02 


4.00C-04 


- 


- 


- 


23.24 


1.58 


^•'J' -0.07 


1.44e-03 


7.00C-05 


1 is+"-ii 
^■-■-8-0.11 


1.10 


860 


49.07 


- 


9 79 + 0.12 
^•'^-0.12 


1.99e-03 


9.00e-05 


- 


- 


- 


89.74 


- 


50+0.29 

ci.D9_Q 26 


2.56e-03 


8.00e-05 




- 


- 


161.40 


- 


4 CQ+O.IS 
*-O'J_0.17 


6.34e-03 


1.40e-04 


- 


- 


- 


6.97 


5.46 




1.21C-03 


1.60C-04 


"■^"-0.04 


1.28 


923 


13.74 


- 


-"-•^"-o.os 


5.94C-04 


8.50C-05 




- 


- 


21.12 


- 


^•y'J-0.12 


9.12C-04 


9.00C-05 


- 


- 


- 


30.48 


- 


9 09+0.17 


1.04e-03 


1.17e-04 


- 


- 


- 


40.64 


- 


9 fic;+0-34 


9.75e-04 


2.39e-04 


^^^^ 


- 


- 


59.76 


- 


9 ot.+0.08 
^•■^■^-0.06 


8.46C-03 


2.80C-04 




- 


- 


15.57 


4.80 


^•'-'^-0.12 


2.64C-03 


3.10C-04 


^'■■'"-0.23 


1.12 


928 


28.11 


- 


9 1:4+0. 28 


3.95C-03 


4.40C-04 


0.69+5['?i 


- 


- 


41.78 


- 


cc+0.45 
o.dd_Q 42 


3.90e-03 


4.90e-04 


1 19+0:51 


- 


- 


57.35 


- 


, (,,+0.43 


5.85e-03 


5.50e-04 


0.761;;:^;: 


- 


- 


78.24 




4 07+0.92 


4.91C-03 


6.40C-04 


o-8ol!!:S 


- 


- 


94.95 


- 


5.i9l»:ll 


1.20C-02 


6.00C-04 


0.48l«:i« 


- 


- 


20.86 


1.17 




1.96C-03 


9.00C-05 


n fiQ+0 09 


1.06 


2144 


33.81 


- 


4 11+0.66 


2.63e-03 


1.30e-04 




- 


- 


46.04 


- 


q 7n+0-55 
•^•'^-0.58 


3.70e-03 


1.60e-04 


- 


- 


- 


59.71 


- 




3.68C-03 


1.70e-04 


51+0-09 


- 


- 


74.82 




4-49lS!;^f, 


4.90C-03 


2.20e-04 








93.52 




6.07+i-t^ 


4.29C-03 


2.00C-04 








114.38 




4.93lj:?« 


4.82e-03 


2.00e-04 


42+0-07 
"-'•*^-0.07 






139.56 






7.44e-03 


2.10e-04 








171.22 




5.79+l-';2 


4.70C-03 


2.30C-04 








215.82 






7.95C-03 


2.40C-04 








9.52 


3.15 


9 CI +0.35 
^•'^^-0.28 


5.70e-04 


4.70e-05 


0.75t°;fo 


1.10 


2756 


15.58 




6 S5+2-48 


9.85e-04 


7.50e-05 








20.78 




, (,-+1.17 

-0.68 


6.48e-04 


7.80e-05 








34.63 




7 qq+1'37 
'•■^^-0.99 


3.24C-03 


1.30C-04 








48.48 




''•'^-'--1.00 


3.19C-03 


1.60C-04 








62.34 




6.96+^-5*^ 


4.22C-03 


2.10C-04 


0.49t°;i| 






77.92 




6.52l-«« 


4.51e-03 


2.20e-04 








95.24 






5.32e-03 


2.90e-04 


0.65t°;?« 






115.15 






3.72e-03 


2.80e-04 
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Table 3 — Continued 



Na.1116 


R t 






Norm 


^Norm 


z 


y2 . 

^rea 


d.o.f 




kpc 


102O cm-2 


kcV 








Zq 






(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


(9) 




138.53 






5.77C-03 


3, 


.20C-04 


U.dU_g 20 








167.10 




9.95tlfr 


6.34C-03 


2 


.80C-04 










202.60 




8.43llf| 


5.79e-03 


4, 


.40e-04 










259.74 




q 02+°-®* 


1.55e-02 


3 


.OOe-04 


(1 7C-+O.O8 






A2052 


13.90 


2.85 


-I 1O+0.10 
^■^^■^-0.08 


3.59C-04 


7, 


.20e-05 


(1 4Q+0.09 


1.21 


1545 




20.22 






1.80C-03 


1, 


.30C-04 










26.96 






1.95C-03 


1, 


.OOc-04 


1 02+°-12 
J^-U^-O.ll 








35.38 




9 09+0.17 


2.28e-03 


1, 


.20e-04 










46.75 




(,1+0.28 


2.04e-03 


1, 


.lOe-04 


(1 74+0.12 
"•'^-0.10 








59.39 




q 1 n + 0.29 

"^■^■^-0.27 


2.43e-03 


1, 


.30C-04 










73.71 




4.oil?,:g 


1.99e-03 


1, 


.OOc-04 


0.47t»:«i 








96.88 




q qq + 0.11 
0.00_Q 


9.19C-03 


1, 


.90C-04 








Hydra A 


19.55 


4.84 


9 Q9+0.19 
^•^^-0.18 


3.17e-03 


1, 


.20e-04 


o.48in^ 


1.16 


999 




33.59 




q (19+0.26 
"^•"^-0.23 


3.87e-03 


1, 


.50e-04 










49.47 




Q oo~f"0.2S 
"^■•-'^-0.26 


5.32C-03 


1, 


.50C-04 


(1 OO+0.04 

U.OO_Q Q3 








72.69 




n 4(5+0.33 


4.93C-03 


1, 


.50C-04 










122.16 




q 7-1+0.14 


1.26C-02 


2, 


.OOe-04 








PKS0745-191 


31.75 


43.39 


04+0.19 


7.42e-03 


2, 


.80e-04 


(1 t;i+0.07 
""^^-0.08 


1.09 


1472 




57.83 




r. 01+0.47 
'J-5J^_0.41 


9.03e-03 


2, 


.80e-04 










91.85 




f, -(,+0.68 


1.07e-02 


3, 


.OOc-04 


(1 or;+0.05 








147.42 




-7 -1 9+0.56 
' ■-'-^-(1,51 


1.38C-02 


3 


,00c-04 










:-! 10.20 




-U,4i 


2,32(--02 




,noc-oi 
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Table 4. Results for 2-T Deprojected, A^h Fixed MEKAL Fits 



Name 






Hot 


Norm 


^ O'PTTh 


TCool 


Norm 






z 




d.o.f 




kpc 


1020 cm-2 


keV 


Hot 


Hot 


keV 


Cool 


Cool 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




(10) 


(11) 


(12) 


2A0335+096 


12.44 


18.11 




1.48e-03 


3.00e-04 


l-06Io.o9 


6.02e-04 


1.09e-03 


0. 


QO+0-20 


1.32 


1968 




19.07 


- 


"'-0.08 


2.87C-03 


3.40C-04 










- 


- 


- 




25.71 


- 


1 Q7+().08 
-0.08 


6.18C-03 


3.80C-04 








0, 


Q-i +0.05 


- 


- 




33.17 


- 


9 o-, +0.12 
^•^-^-0.11 


5.27C-03 


2.90C-04 










- 


- 


- 




42.29 


- 


O yii+0. 15 
^- '•J-0.14 


6.34e-03 


2.70e-04 








0. 


Sfi+0-04 
■°"-0.07 


- 


- 




54.31 


- 


2 q5+"'^^ 


6.50e-03 


2.50e-04 










- 


- 


- 




70.48 


- 


3.62t»;^?, 


5.76C-03 


2.30e-04 










- 


- 


- 




90.80 


- 




8.58C-03 


2.30C-04 








0, 


f;c: + 0.05 
■"'^-0.05 


- 


- 




165.84 




A zL^+0-11 


2.50C-02 


4.00C-04 


















D.y ( 


u.4tD 


, fi4+0.22 




y. / ue-uo 


on+O-03 




i.oUe-U4t 


1 

1 . 


C9+I.34 
■°^-0.57 




001 




13.74 




-'-•^■^-0.04 


5.80e-04 


7.30e-05 








1. 


PC -1+0.13 
■°^-0.18 








21.12 




1 90+005 

J^-y'-'-o.ii 


9.42C-04 


7.10C-05 


















30.48 




^•'^-'--0.13 
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Table 5. Fits to Temperature Profiles from Deprojected X-ray Spectra 



Name N Tioo Slope Prob 

(keV) 
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Table 6. High Resolution Entropy Fit Results (Evenly Binned 5 Arcsecond Bins) 
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Method 


(keV cm^) 


(keV cm^) 


Slope 




Prob 


zAUooo 


QQ 
OO 


U.io 


i 


o.4oo ± U.o4 


1 A _i_ O 

ioo.4 ± z.o 


1 /1 1 -1- n noi^ 
i.4i ± U.UzD 


i^i 
bi 


1 /-, 

ie-o 








i 


= U 


1 on ^ 1 o o 
izy.D ± z.o 


1 1 /I -1- n ni 
i.i4 ± U.UiO 


Zio 


n 
U 








o 
z 


/.z4U ± U.oo 


■to A 7 _L O c; 
io4. / ± z.o 


1 /I n -L n noQ 
i.4y ± U.Uzo 


1 A 

/4 


0^ nc; 

oe-Uo 








o 
z 


n 

= U 


1 Ofi -| 1 o o 
iZO. i ± Z.O 


i. iU ± U.Uio 


4bo 


U 


A 1 QQ 

Aioo 


zo 


U.io 


i 


1 1 QO 1 1 1 

ii.oz ± i.i 


1 /I n 1 1 q Q 

i4y. i ± o.o 


1 Qi -J- n nc^n 
i.oi ± U.UoU 


Z I 


n nn 

u.uy 








i 


= U 


1 t^n 1 1 o r 
iDU.i ± O.D 


U.yb 1 ± U.Uzi 


1 1 o 

iiz 


O^ 1 /I 

ze-i4 








o 
z 


io. /U ± i.U 


14o.i ± o.o 


1 on _i_ r\ r\K A 
i.oy ± U.Uo4 


OO 
OO 


A AO 

U.Uz 








o 
z 


A 

= U 


1/1/1 O _1_ Q 

i44.o ± o.o 


n nn^^ -L n non 
U.yUb ± U.UzU 


ioo 


n 
U 


A 1 7nc; 
Ai i\)o 


Q1 


U.io 


i 


1 A K'7 1 -I T 

i4.o/ ± i./ 


1 on o 1 o o 
izU.o ± o.o 


1 1 Q -L n n<50 
i.io ± U.Ubo 


O.o 


n nn 

u.yy 








i 


= U 


1 OO 0-1-0 

ioz.o ± o.o 


n Q /i o -L n noc: 
U.o4z ± U.Uzo 


c;i 
Oi 


n nno 
U.UUo 








o 
z 


iy.y/ ± i.D 


11/1 O -L Q 

ii4.z ± o.o 


1 on -L n n^;n 
i.zy ± U.Uby 


n 

y.o 


n nn 

U.yy 








o 
z 


A 

= u 


1 on o -L 
izy.z ± o.o 


n 7Q/I -L n no/i 
U. / o4 ± U.Uz4 


1 nn 
iUU 


0/~. 1 n 
ze-iU 


A onon 
AzUzy 


OQ 
ZO 


n on 
U.zU 


i 


D.o4/ ± O.O 


1 ^O /I _1_ o 

ioz.4 ih O.z 


n Qni _u n n^^Q 
U.oyi ± U.Ubo 


n 
O.U 


n nn 

U.yy 








i 


A 

= u 


loo.o ± O.U 


n vno 1 A AO /I 
U. / yo ± U.Uo4 


/ .0 


n nn 

U.yy 








o 
z 


iU.io ± d.o 


lOO.l ± O.Z 


n noQ _i_ n nvn 
U.yzo ± U.U lU 


o.y 


n nn 

U.yy 








o 
z 


A 

= U 


lb / .0 ± O.U 


n /I _i_ n noo 
U. / b4 ± U.Uoo 


iU 


n nn 

U.yy 


AzUOz 


on 
zU 


U.U iU 


i 


n 001 _L n CO 
y.zzi ± U.oz 


1 A 1 A ^ 

ioo.u ± y. / 


1 _l n n^^ 
i.ob ± U.Ubb 


7n 
/U 


o^ nn 

ze-uy 








i 


A 

= u 


1 KO c: -L c; n 

ioo.o ± o.y 


1 1 n -L n non 
i.iU ± U.Uzy 


i / b 


n 
U 








o 
z 


1 1 TO _i_ n TQ 

ii. 1 6 ±0.10 


1 nT A 1 1 1 

iy / .y ± ii. 


1 Tn -1- n nTc; 
i./U ± U.U/0 


O/ 


1^1/1 
ie-i4 








o 
z 


A 

= U 


i4 / .y ± o.o 


1 n7 -1- n non 
i.U/ ± U.Uzy 


O70 
Z / O 


n 
U 


A262 


31 


0.050 


1 


1.022 ±0.38 


230.0 ±8.8 


1.11 ±0.028 


154 











1 


= 


218.2 ±6.6 


1.05 ±0.013 


162 











2 


3.606 ± 0.34 


241.6 ± 10. 


1.20 ±0.031 


190 











2 


= 


195.9 ±5.9 


0.979 ±0.013 


290 





A496 


26 


0.080 


1 


7.109 ± 1.0 


175.0 ±9.6 


1.20 ±0.065 


5.8 


0.99 








1 


= 


150.8 ±6.6 


0.902 ± 0.029 


39 


0.01 








2 


11.32 ±0.97 


177.2 ± 10. 


1.31 ± 0.073 


8.2 


0.99 








2 


= 


133.0 ±5.8 


0.788 ± 0.027 


90 


3e-10 


Hydra A 


20 


0.10 


1 


13.43 ± 1.1 


98.74 ± 3.9 


1.22 ± 0.078 


2.2 


0.99 



-44- 



Table 6 — Continued 



Name N Radius Method Kq -f^ioo Slope Prob 

(Mpc) (keV cm^) (keV cm^) 
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Table 7. High Resolution Entropy Fit Results (Evenly Binned 2.5" Annuli) 



Name 
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Radius 
(Mpc) 


Method 


(keV cm^) 
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Table 7 — Continued 



Name N 


Radius 
(Mpc) 


Method 


Ko 
(keV cm^) 
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Slope 
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Table 8. Entropy Fit Results (Evenly Binned 2.5" Annuli, Temperature from best-fit 

deprojected temp) 
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